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1. INTRODUCTION 

In continuum and fluid dynamical models, particles, which leave the system and reach 
the detectors, can be taken into account via freeze-out (FO) or final break-up schemes, 
where the frozen out particles are formed on a 3-dimensional hypersurface in space-time. 
Such FO descriptions are important ingredients of evaluations of two-particle correlation 
data, transverse-, longitudinal-, radial- and cylindrical- flow analyses, transverse momen- 
tum and transverse mass spectra and many other observables. The FO on a hypersurface 
is a discontinuity, where the pre FO equilibrated and interacting matter abruptly changes 
to non-interacting particles, showing an ideal gas type of behavior. 
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The frequently used Cooper- Frye formula, the rapidity distribution, transverse momen- 
tum spectral, etc., - all include expression fpoi^, P, T, n, u^) p'^da^ , where fpoi^, P, T, n, u^) 
is the post FO distribution which is unknown from the fluid dynamical model, da^^ is nor- 
mal vector to the FO hypersurface. Those formulas work well for timelike da'^ [p'^da^ > 0). 
If da^ is spacelike we count particles going backwards through FO front as well as out- 
wards. The post FO distribution can not be a thermal one! In fact fpo should contain 
only particles which cross the FO front outwards, p^da^ > 0, since the rescattering and 
back scattering are not allowed any more in the post FO side. 

fFo{x,p,T,n,u'',daf,) = 0, p^rfa^ < . (1) 

Initially if we know the pre FO baryon current and energy-momentum tensor, Nq and 
Tq'^ , we can calculate locally, across a surface element of normal vector (ia^ the post 
FO quantities, A^^ and Ti"", from the relations: [A^'' da^,] = and [T'"' da^] = 0, where 
[A] = A — Aq. One should also check that entropy is nondecreasing in FO. In numerical 
calculations the local FO surface can be determined most accurately via self-consistent 
iteration. Initial ideas to improve the Cooper-Fry FO description this way were suggested 
inrefs. [HHI. 

2. FREEZE-OUT DISTRIBUTION FROM KINETIC THEORY 

Let us assume an infinitely long tube with its left half (x < 0) filled with nuclear matter 
and in the right vacuum is maintained. We can remove the dividing wall at t = 0, and then 
the matter will expand into the vacuum. Continuously removing particles at the right end 
of the tube and supplying particles on the left end we can establish a stationary flow in 
the tube, where the particles will gradually freeze out in an exponential rarefraction wave 
propagating to the left in the matter. We can move with this front, so that we describe 
it from the reference frame of the front (RFF). 

We assume that there are two components of our momentum distribution: ffree{x,p) 
and fint{x,p)- However, we assume that at x = 0, ffree vanishes exactly and ftnt is an 
ideal Jiittner distribution, then fint gradually disappears and ffree gradually builds up, 
according to the differential equations: 

d^fintix,p)dx = -Q{p''da^) ^"^^^ /mt (x, p)dx , 

(2) 

dxffree{x,p)dx = +&{p^ dci f,) ^^^^ ^ fint{x , p) dx , 

where cos^p- = ^ in RFF. It expresses the fact that particles with momenta orthogonal 
to the FO surface leave the system with bigger probability than particles emitted at an 
angle. 

Such a dramatically oversimplified model can reproduce cut Jiittner distribution as 
a post FO one, but it does not allow the complete FO - the interacting component of 
momentum distribution survives even if x — > oo. 

To improve our model we take into account rescattering within the interacting com- 
ponent, which will lead to re-thermalization and re-equilibration of this component. We 
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use the relaxation time approximation to simplify the description of the dynamics. Thus, 
the two components of the momentum distribution develop according to the differential 
equations: 

dxfintix,p)dx = -Q{p''da^)^^fintix,p)dx+ 

(3) 

+ [feq{x,p) - fintix,p)] jrdx, 
d^ffree{x,p)dx = +Q{p'' dci ^Y^^ fint{,X , p)dx . (4) 

The interacting component of the momentum distribution, described by eq. (^, shows 
the tendency to approach an equilibrated distribution with a relaxation length A'. Of 
course, due to the energy, momentum and particle drain, this distribution feg{x,p) is not 
the same as the initial Jiittner distribution, but its parameters, neg{x), Teq{x) and u^g{x), 
change as required by the conservation laws. 

Let us assume that A' ^ A, i.e. re-thermalization is much faster than particles are 
freezing out, or much faster than parameters, neq{x), Teq{x) and u'^q{x) change. Then 
fint^x.'p) ~ /eg(x,p), for A' -C A . For feq{x,'p) we assume the spherical Jiittner form at any 
X including both positive and negative momentum parts with parameters n{x), T{x) and 
u^pq{x). (Here Up^pQ^x) is the actual flow velocity of the interacting, Jiittner component, 
i.e. the velocity of the Rest Frame of the Gas (RFG) In this case the changes of 
conserved quantities due to particle drain or transfer can be evaluated for an infinitesimal 
dx and new parameters of feq{x + dxjp) can be found 0. 

We would hke to point out that, although for the spherical Jiittner distribution the 
Landau and Eckart flow velocities are the same, the change of this flow velocity calculated 
from the loss of baryon current and from the loss of energy current appear to be different 
duf^p^jiPQ^x) ^ du^ p j^pQ^x) . This is a clear consequence of the asymmetry caused by the 
FO process as it was discussed in ref. 0J^. This problem does not occur for the freeze-out 
of baryonfree plasma. 

We performed the calculations, according to this model, for the baryonfree and massless 
gas 1^,0]. We would like to note that now fint{x,p) does not tend to the cut Jiittner 
distribution in the limit x —>■ oo. Furthermore, we obtain that T — * 0, when x — > oo. So, 
fin.t{x,p) = -^2^^^ exp[(/i — p'^Uu)/T] — > 0, when x ^ oo. Thus, all particles freeze out 
in the improved model, but such a physical FO requires infinite distance (or time). This 
second problem may also be removed by using volume emission model discussed in [0. 
The application of our onedimensional model to transverse expansion gives result showing 
in Fig. 1, which is in qualitative agreement with experiment. 

Recent calculations have confirmed that the FO hypersurface idealization can be jus- 
tified even in microscopic reaction models (like UrQMD or QGSM) for nucleon data in 
collisions of massive heavy ions. 

The improvements presented here are essential and lead to non-negligible qualitative 
and quantitative changes in calculations including FO. Several further details and conse- 
quences of this improved approach have to be worked out still (e.g. |^) to obtain more 
accurate data from the numerous continuum and fiuid dynamical models used for the 
description of heavy ion reactions. 
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Figure 1. The local transverse momentum (here px) distribution for baryonfree, massless 
gas at Py = 0, X = lOOA and Tq = 130 MeV. The transverse momentum spectrum is 
obviously curved due to the freeze-out process, particularly for large initial flow velocities. 
The apparent slope parameter increases with increasing transverse momentum. This 
behavior agrees with observed pion transverse mass spectra at SPS [§]. From P]. 
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